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Abstract 

We show how singularities shape the evolution of rational discrete dynamical sys¬ 
tems. The stabilisation of the form of the iterates suggests a description providing 
among other things generalised Hirota form, exact evaluation of the algebraic entropy as 
well as remarkable polynomial factorisation properties. We illustrate the phenomenon 
explicitly with examples covering a wide range of models. 


1 Introduction 


Many of the algebraic aspects of rational discrete systems have already been investigated, 
especially in view of their integrability. A nnmber of these aspects are extensions to the 
discrete case of featnres of continnons systems, bnt the rationality of the evolntion made 
the algebro-geometric approach inescapable. It is indeed at the basis of any classihcation 
attempt, symmetry or mnltidimensional consistency analysis, as well as complexity measnre 
via algebraic entropy. See for example the nnmerons results exposed in the series of SIDE 
meetings (^. 

The nature of the evolution is also responsible for one recurrent fact: looking at a h- 
nite number of steps of the discrete evolution yields informations which are in essence of 
asymptotic nature, like integrability, hierarchies, or value of the entropy. 

Motivated by the original works and comforted by the more modern approaches 
^ on continuous systems, the importance of the singularity structure was recognised very 
early [^. The use of the apparatus available in two dimensions, and notably the theory of 
intersection of curves on algebraic surfaces, then lead to powerful theorems, in particular on 
discrete Painleve equations and “QRT” maps [8 10 


Direct computations of discrete evolutions have also been performed, especially to detect 
integrability, endeavouring to reduce the size of the calculations. For instance, looking at 
the images of a straight line in the space of initial conditions, inspired by the geometrical 


idea of 11 , allows to produce an exact sequence of degrees of the iterates, and in turn to 
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12,13 


Restricting the evolution to integers lightens 

gives an 
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evaluate exactly the algebraic entropy 

even more the calculations. Looking then at the growth of their height 
approximate but efficient way to evaluate the algebraic entropy. Going even further one may 
perform the calculations on hnite helds, analyse various statistical properties 16,17 , and 
eventually detect integrability. 

We take here an opposite attitude, and choose an ingenuous option: we do evaluate 
exactly the hrst steps of the evolutions, and analyse their structure, and especially their 
factorisation properties. 


The main outcome is that the form of the iterates, shaped by the singularities, suggests 
changes of description which automatically provide: 


Exact calculations of the algebraic entropy, and proofs of its algebraicity (see 18 20 
for a general point of view on this question) 


Generalisations of the discrete Hirota-Sato form and r functions 21-27 , reducing to 
the standard quadratic form in specihc integrable cases 

Various polynomial and integer factorisation properties similar the Laurent prop¬ 
erty 


28 33 


The plan of the paper is the following: 

In section we briefly recall basic algebro-geometric notions which will be at the core of 
the phenomenon we exhibit. 

In section we describe explicitly a number of models in various dimensions: 

1 - An algebraically integrable two dimensional map in the QRT family: McMillan 

2 - A discrete Painleve equation dPn, i.e. a non autonomous extension of the previous 

3 - Another discrete Painleve equation: qPyi 

4 - A non integrable non conhning map in two dimensions: Jaeger 

5 - A conhning non integrable map in two dimensions: JNH-GMV 

6 - Another conhning non integrable map in two dimensions 

7 - A three dimensional algebraically integrable map: N=3 periodic Volterra 

8 - A linearisable map 

9 - An unruly model in three dimensions 

10 - A recurrence of order 2 on functional space. Delay Diherential equation 

11 - An integrable lattice map: Qv 


We conclude with suggestions for further explorations. 


2 



















2 The importance of being singular 

We use complex projective spaces as spaces of initial conditions X. Suppose for simplicity 
that X is of dimension N, with iV + 1 homogeneous coordinates and call ip is the forward 
map, and ip the backward map. Then 

p : [xo, xi,..., Xn] -t [xq, x\, xjv] 

^ : [Vo, 2/1, • • •, Vn] -t [Vo, y[, • • •, 2/^] 

The evolution step will always be given by a birational map, so that both the forward and the 
backward evolution are described by polynomial maps, i.e. x'-, j = 0 ... n and i/', j = 0 ... n 
are polynomials of degree and dp respectively (usually = dp > 1). 

The singular points of p (resp. ip) are the ones for which = 0,j = 0...n, (resp. 
y'j = 0,j = 0...n). We know that the sets of singular points are algebraic varieties of 
dimension < N — 2. 

Since ‘ip is the inverse of p' means that the composition p ■ ip appears as a multiplication 
of all coordinates by a common factor, we have the two basic relations 

Ip ■ P ■ id, P ■ Ip ■ id (1) 

The two polynomials and Kp, both of degree d^dp — 1, may be decomposable. 

= = ( 2 ) 

i=i j=i 

Each factor dehnes an algebraic variety of codimension 1 playing an important role in 
the sequel. 

There is a simple relation between the varieties and Kp and the singular locus of p 
and ip\ the varieties of equation = 0 are blown down by p and their images are entirely 
made of singular points of ip. This reflects the fact that one cannot take one step forward 
and then a step backward when starting from a point on = 0. The same applies to ip 
mutatis mutandis. 

Suppose S is an indecomposable variety of codimension 1 of equation E^, = 0. The 
pullback by p of the equation of S gives the equation Ej^/ of the image S' of S by '0- The 
important point is that this pullback may contain additional factors. 

P*{Ej:) = E^, 

Such factors are necessarily built from components of Their presence reflects the differ¬ 
ence between total and proper transform: the non-singular points of S go into the proper 
transform. The singular subvarieties contained in S are blown up to the other components. 
In a way the proper transform is the true image, disregarding the singularities. 

Singularity conhnement in this context is just that some iterate of p sends components 
of Kip into components of Kp. The regularisation comes, when described with homogeneous 
coordinates, from the removal of factor common to all coordinates. This is also the origin 
of the possible drop of degree of the iterates of p. 

Remark 1: A consequence of the existence of the variety = 0 is that the image by p 
of a generic line always hits some singular point of ip. The reason is that any generic line 
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crosses since we are working with complex projective space. In other words, while the 
image of a generic point is always a generic point, the image of a generic line is never a 
generic line. 

Remark 2: A variety could be its own transform. It is then covariant by (p and 'ip. 


Covariant objects play a fundamental role in the description of the algebraic invariants 34 


By abuse of language we will say that S' is the proper transform of S', if S and S' are 
the equations of varieties which are the proper transforms of each other by (p or 

Given an evolution map (p, denote by pk the successive image^of po = [xq, xi,..., x^]- 
The components of pk will factor into indecomposable blocks. These blocks are necessarily 
either the factors of the multiplier K^p and their transforms, either the transforms of the 
coordinate planes. These blocks verify remarkable algebraic recurrence relations, and this is 
the subject of this paper. 

Remark 3: Birational changes of coordinates, which dehne the natural equivalence rela¬ 
tion between different descriptions of the same model, affect the singularity structure, the 
value of Ky, and and the form of the equations relating the various blocks. They however 
will not spoil the general features of the recurrences between blocks. 


3 Eleven models 

The simplest possible type of discrete systems is given by recurrences of hnite order. A 
recurrence of order k may be looked at as a map in its k dimensional space of initial data. 
There are two natural generalisations, leading to inhnite dimensional space of initial condi¬ 
tions: recurrences dehned over functional space, and recurrences with multi-indices (lattice 
maps). Both will be considered, again supposing rational invertibility of the evolutions. 

This section contains the explicit description of the aforementioned blocks and recurrence 
relations for eleven different models, integrable as well as not integrable, hnite dimensional as 
well as inhnite dimensional, to oher a panoramic view on the property we describe, including 
a limiting case (section [T^ . 

3.1 McMillan 

The model is a prototype of algebraically integrable map in two dimension, belonging to the 
Quispel-Roberts-Thompson family [8 

The map ip associated to the model reads 

(p : [x, y, z] — > [—y (x^ — z^) + 2 axz"^, x (x^ — z (x^ — (3) 

Its inverse is 

ip : [x, y, z] — > ly(y- z) (y + z), xz"^ -y‘^x + 2 yaz"^, z{y- z){y + z)] (4) 

and 

f) ■ (p = {x - zf {z ^ x)^ = p ■ f) = {y - zf {y + zf (5) 

^The point Pk+i can be obtained from pk either by the action of (/? on pk, either by pulling back the 
coordinates of pk. The homogeneous coordinates obtained in these two ways may differ, but they represent 
the same point projectively. 
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The form of the hrst iterates 


Po = [x,y,z\ 

Pi = [Ai, xBiCi, zBiCi] 

P 2 — [^2-BlCi, AiB2C2-i ^-Bi C*! i?2 (^ 2 ] 

( 6 ) 

Pk = [Ak Bk_iCk-i, Ak-i BkCk, z Bk-iCk-i BkCk] (7) 

Expressing that pk+i = p{pk) gives only one condition: 

{Ak — zBkCk) {Ak + zBkCk) {Ak-iBk+iCk+i + Ak+iBk-iCk-i) 

— 2 a z"^ Ak Bk-i Ck-i Bk Ck Bk+i Ck+i = 0 . ( 8 ) 


This condition does not suffice to determine {Ak+i, Bk+i,Ck+i}. 

Claim: There exist algebraic relations between the blocks A, B,C. These relations allow 
to calculate {Afe+i, -Bfc+i, Ck+i} in terms of the previous A, B, Cs . Moreover A^+i, Bk+i, Ck+i 
are the proper transforms of Ak,Bk,Ck- 
Proof. We have, for k > 3: 


Ak — z BkCk + Bk+iCfc-i — 0 , 

^fc-iBk+iCk+i + Ak+ii?fc-iCfc_i 


Ak + z BkCk — 7?fc_iCk+i — 0 
V 2 a z^ Ak Bk Ck = 0 


(9) 


Equations (|^ imply ([^, and can be verified directly for k=3 and k=4. The validity for 
general k is obtained by recursion. Proving that the form of Q and of relations ([^ is stable 
is just a matter of counting factors i?’s and C’s. We know from section ([^ that 

p*{Ak) = Ak+i 

p*iBk) = Bk+i 

p*{Ck) = Ck+i 

for some exponents a, /?, 7 . Using ([^ we get Ak+i from the previous A, B, C’s. We may thus 
evaluate all the exponents a{k + l),(3{k + l),'y{k + 1). The outcome is that Ak+i, Ak+i, Bk+i 
are the proper transforms of Ak, Ak, Bk and do not factorise, as no new factors Bi or Ci are 
left over in the components of pk after k = 3. QED. 

We will present similar properties in the subsequent sections. Their proof goes along the 
same lines and will not be detailed. 

Relations ([^ define completely the evolution oi {Ak,Bk,Ck}- They extend over a string of 
points of length 3. Although their solution is written as fractions, the result is automatically 
a polynomial in terms of the initial eonditions [x,y,z]. They moreover enjoy a Laurent 
property 


28 31 33 


Dehne a map ^ : [P, Q, R, U, U, W] —^ [U, U, W, U', V, W'] with 

-PV'W'-2az‘^UVW U-zVW U + zVW 

U = -——-, V = ---, W = - 


QR 


R 


Q 


^There always is a rescaling possibility of the various factors, in particular an ambiguity in the signs, 
together with a possibility of exchanging B and C. 
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which implements the solution of ([^ as a map. We may consider iterations of starting from 
arbitrary initial data [p, q, r, u, v, w]. The images are Laurent polynomials in [p, q, r, u, v, w]. 
If in addition the triplets \p,q,r] and [u,v,w] happen to be of the form [Ak-i,Bk_i,Ck-i\ 
and [Ak, Bk,Ck], then the iterates are polynomials. 

Setting Tfc = BkCk, one may rewrite the ’raw’ form ([^ as 

I Al-z^Tl + Vk-iVk+i = Q 

I ^fc-iTfc+i + ^fc+iTfe-i + 2 a = 0 


which is nothing but the bilinear quadratic discrete Hirota form, the inhomogeneous co¬ 
ordinates of the kth iterate being just [Ak/{zVk),Ak-i/{zVk-i)], as suggested in 

following (^[^. 

The recurrence Q also gives the constraints obeyed by the sequence of degrees (5(A„),... 
of the successive A, B,CA, 
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I ^(^fc+i) + ^{Bk-i Ck-i) — S{Ak-i) + 6{Bk+i Ck+i) — S{Ak) + S{Bk Ck) + 2 

[ ^{Bk+i) + = ^(C'fc+i) + ^{Bk-i) = 5{Ak) = S{Bk Ck) + 1 

and consequently the one verihed by the degree dn = d{An) + 6{Bn-i Cn-i) of p„ 


dn ^ d^—\ “ 1 “ 3 dn—2 dyi—^ 0 , 


which proves quadratic growth of {dn} and vanishing of the algebraic entropy. 
Finally the invariant of the model may be rewritten 


I 


A ^ _j_ / Ak+i \ ^ / ^kAk+i _ ^ 

zTk) \zTk+iJ \z‘^TkTk+i ) 


3.2 dPn 

This an integrable non autonomous extension of the previous model, and has the Painleve II 
equation as a continuous limit 39 . Its non autonomous nature invites us to write the map 
in three dimensions, the added variable having a linear evolution. 

if : [x, p, z, t) —)■ [—y{ld — x^) + {cz + ht) xt,x {C — x^) ,{z + t) {C — x^) , t — x^)] 

and 

V’ : [a;, P, t] —^ [p - p 2 ) , -x{t^ -y‘^)+t {bt + c(z - t)) p, {t - z) {C - p^) , t - p^)] 


= {t + x)'^ {t - x)^ = B{ = (t -f y)^ {t - yf 


The sequence of point we get from po = [x, p, z, t] is 


Pi = [Ai, x^iCi, (z + t) BiCi, tBiCi] 

P2 = [A2 Bi Cl, Ai B2 C2, {z 2 t) Bi Cl B2 C2, t Bi Cl B2 C2] 


Pk — [Ak Bk-iCk-i, Ak-i BkCk, {z + kt) Bk-iCk-iBkCk, t Bk-iCk-i BkCk] (10) 
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The recurrence on Ak Bk Ck is given by the following constraints, generalising straightfor¬ 
wardly ([^; 


Ak — tBkCk + Bk+i Ck-i — 0 , Ak + tBkCk — -Bfc-iCk+i — 0 

Ak-i Bk+i Ck+i + Ak+i Bk-i Ck-i — t (bt + c (z + kt)) Ak Bk Ck = 0 


( 11 ) 


The proof is similar to the one given in the previous section. Setting here again Tk = BkCk 


we get the ‘Hirota form’ found in 37 , extending over a string of length 3. 


We can use these relations to prove vanishing of the entropy, and check that the vari¬ 
ous Ak,Bk,Ck are the proper transforms of Ak-i, Bk-i,Ck-i. The algebraic invariant has 
disappeared, but the overall algebraic structure is essentially unchanged, compared to the 
previous model, apart from one coefficient which became non-constant. 


3.3 qPvi 

The map 93 , which is a discrete version of the Painleve equation Pyi, may be written as the 
composition three maps, taken from equations (19,20,21) of |40] . 

^(Po) = y, Z, f]) = (P 3 • (^2 • y, Z, f]) with 

Pi{Po) = [xy [x — ct) [x — dt) ,sh{x — az) [x — hz) C, zy {x — ct) [x — dt) ,ty {x — ct) [x — dt)] 
P2iPo) = [cd {y - pz) {y - zr) C, yx {-st + y) {y - ht ), 2:0; {-st + y) {y - ht ), tx {-st + y){y- ht)] 

PsiPo) = [x,y,qz,t], q= {cdpr)/{absh) 

The hrst iterates yield the following sequence of points, setting a = ab, (3 = cd, a = sh, 
and p = pr: 

Pi = [a (3 y Ai Bi Ei Fi, a P xC Ci Di Gi Hi, fdpxyz Ci Di Ei Fi, aaxyt Ci Di Ei Fi] 

P2 = xt^ A2 B2 E2 F2, ci /3 P y C2 D2 G2 H2, ildp)"^ z C2 D2 E2 F2, (crcT t C2 D2 E2 T2] 

P3 = [(a/5) y As B-i E3 F3, a /5 x C3 D3 G3 H3, {jdpf xyz C3 D3 E3 F3, {aaf xyt C3 D3 E3 F3] 
P4 = [{a( 3 YP X A4 i?4 i?4 F4, a yt^ C4 ZI4 G4 H^, {IdpY z C/i E^ F4, (aa)^ t C4 ZI4 E^ F4] 

P5 = [(a/ 5 ) yC As B3 E3 F5, a x C5 D3 G5 H3, {jdpf xyz C3 D3 E3 F5, («a)'' xyt C3 D3 E3 F5] 
P6 = [(a/ 5 ) a; Ag B^ E^ F^, a y C^ Gg H^, {^pf z Gg Ee F^, {aaf t Gg Du E^ Fg] 


The form of pk is thus 

Pk = [(«/5) V'’- AkBkEkFk, a f3’^-^Pf2,k CkDkCkHk, (/5p) CkDkEkFk, {aa)’^U,k CkDkEkFk]{l2) 

where the fi^k depend on the initial conditions, are such that ft^k = fi,k+6, and can be read 
from the iterates given above. They are summarised in the table 


k mod 6 

i = 1 

i = 2 

i = 3 

/ = 4 

0 

X 

y 

z 

t 

1 

y 

xC 

xyz 

xyt 

2 

xC 

y 

z 

t 

3 

y 

X 

xyz 

xyt 

4 

X 

yC 

z 

t 

5 

yD 

X 

xyz 

xyt 
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The various factors A, B, C, D, E, F, G, H verify simple recurrence relations, which can 
be checked for the hrst few ones, and then proved by recursion, as in the previous sections. 
They read 


P ^ ^i,k sh GkFIA.\^Bi^—\ 

P ^ ^i,k Ei^Fj~ shoj\ i^ Gj^Hj; 

1 / 2 ,A: EkFk - U2,k GkHu = C^Dk 

U2,k EkEk - (cd)"-^ U2,k GkHu = Gu-iD^ 

giving Ak,Bk,Gk, in terms of A^-i, Bk-i, Ek,Ek,Gk,Hk, and 


^ Xi^kAkBk — sh fii^kGkDk — Ek+iTfc 
Xi^kAkBk — sh fii kGkDk = EfcFk+i 
a%^-\sh)^-^ X2,kAkBk - (pr)^ P2,kGkDk = 
a’‘-^b\sh)^-^ X2,kAkBk - (pr)^ P2,kGkDk = G^Hk+i 


giving Efc+i, -ffc+i, Gfc+i, idfc+i in terms of Ak, B^, Gk, D^. The coefficients A, /r, i /,oj appearing 
in the previous relations depend on the initial conditions and are given by the following table: 


k mod 6 

vi,k 

1^2,k 

^l,k 

^ 2 ,k 

Ai,fc 

X 2 ,k 


p 2 ,k 

0 

z 

t 

y 

y 

X 

X 

t 

Z 

1 

yz 

y 


t 

1 

1 

t 

xz 

2 

z 

t 

y 

y 

xt 

xt^ 

1 

z 

3 

y 

yt 

1 

1 

1 

1 

xt 

xz 

4 

z 

1 

yP 

yt 

X 

X 

t 

z 

5 

y 

yt 

1 

1 

t 


X 

xz 


As in the previous case, equations ( l4p!3 ) can be put in a quadratic form by setting 

Uk — Af^ B]^^ Vfc ~ Cfe hhfc = Ef^ Fjf^ Tj« = G]^ 

with coefficients depending on the initial conditions and on k. They differ from the ones 
found in 
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Relations (13) and (14) define an iteration, starting from [Ai, Ri, Ci, Hi, Ei, Fi, Gi, idi]. 
The explicit calculation of the first few iterates indicates that they are Laurent polynomials 
in the initial conditions. If in addition the initial conditions were produced by the action of 
the map ip on some po = then all resulting quantities are polynomials in x, y, z, t. 


Relations (13) and (14), as well as (12) may also be used to prove vanishing of the entropy. 


because they provide an exact evaluation of the successive degrees of the iterates of p. 


3.4 A non integrable non confining map in P 2 

The model was proposed by F. Jaeger in relation to studies of Bose-Meisner algebras (see 
for example |^). It is the product of two birational involutions of P 2 -, constructed from a 



























projective linear map L defined by the matrix 

’ 1 2 2 ' 

13- 4 

14- 5 

and the fundamental involution 

j : [x,y,z] —^ [yz,zz,xy]. 


Define 


i = L ^ ■ j ■ L, = i ■ j, xjj = j ■ i, Kip = = 2bxy z. 

It will appear that the pattern leads naturally to decompose ip as the product 

^ = {L-^ j) ■ (L j) 

The sequence of iterates po, qi = = Lj{qi ),... we obtain is 

Po = [x, y, z] 
qi = [Ui,Vi,Wi] 

Pi = [Ai, xBi, xCi] 
q2 = [U2, Ui 14, Ui W 2 ] 

P2 = [A2, Ai B2, Ai ,C2\ 

Qk = [Uk, Uk-iVk, Uk-iWk] 

Pk = [Ak, Ak-iBk, Ak-iCk] 

The recurrence on the blocks A, B,C, U, V, W take now the form 
Ak-i BkCk + 2Ck Ak + 2 Ak Bk = 5 Uk+i 

Ak-i BkCk + ^Ck Ak — 4: Ak Bk = 5 f/fc V^+i 

Ak-i BkCk + 4 :Ck Ak — 5 Ak Bk = 5 Uk Wk+i 


(15) 


(16) 


Uk Vk+i Wk+i + 18 114+1 Uk+i — 14f/fc+i 14+1 = Ak+i, 
Uk Vk+i 114+1 — 7 ITfc+i f/fc+i + 6 Uk+i 14+1 = Ak Bk+i, 

Uk 14+1114+1 — 2 Wk+i Uk+i + Uk+i 14+1 = Afc Ck+i 


(17) 


The recurrence relations we get here are not quadratic. The main point is that in (16) 
and Uk+i and Ak+i have polynomial expressions. This reflects the singularity structure 
of the map: in the iteration process, the lines {y = 0} and {z = 0} are blown down to points 
which whose images never meets any singularity, while the line {x = 0} goes to [1,1,1], 
and then [1,0,0] which is singular. Relations (16) and © allow us to calculate exactly the 
entropy e = log((3 + \/b)/2) of (p. 
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3.5 A confining non integrable map in P 2 

We briefly mention here what has become a prototype of conhning but chaotic map in two 


dimension, described in 42 


if : [x, y, z] —y + az^ — yx‘^, x^, x'^z] 

coming from the simple order 2 recurrence 


^n+l 4“ ^n —1 4“ 9 


Ut 


Here 


Kj(p — X 


This map has been shown to have positive entropy by various methods, among which the 


construction of a rational surface over P 2 where the singularities are resolved 43 . The lift 


of the map to the Picard group of this variety is a linear map whose maximal eigenvalues 
gives the entropy. 

The generic iterate has the form 


Pk — [-^k-3 ^k, ^k-A ^ ^fc-3 ^k-2 
and the recurrence relation between the blocks A becomes 


43 /l3 I ^ ^3 46 a6 _ 43 A a2 _ a2 43 a 
^k^k-3'^^ ^k-l^k-2 ^k-1 ^k-A ^k — ^k-3 ^k-2 


k+1 


(18) 


(19) 


which is equ (4.6) of 30 , extending over a string of length 6 . This relation is not quadratic 
nor multilinear, but it allows to prove (18) with the same type of argument as in the previous 
sections, providing the recurrence condition on the degrees of the iterates of <p, and the value 
of the entropy e = (3 + a/5)/2. 


3.6 Another confining non integrable map in P 2 


Another interesting conhning non integrable map described in 13 , eq 29, was also examined 
for a = 0 in 30 . The map comes from the order 2 recurrence 


Un+l ■ — Un-\ -h 0, 

^n. 


if : [x, y, z] —)■ [z (a;^ 4 - 2 :^ 4 - axz) , yx'^, xyz], 

: [x, y, z] y^x, z 4 - 2 ;^ 4 - azy) , xyz\. 

= x^y'^z {x -\- z{a + a;)/ 2 ) [x + z{a — a)/ 2 ) = x^y‘^z Bi Ci, 

= x^y^z {y + z {a + a)/2) {y + z {a — a)/2) 

with a = — A. 

By blowing up 18 points, we may dehne a rational surface over P 2 where the lift of the 
map becomes a diffeomorphism, and the entropy is given as the logarithm of the inverse of 
the root of — 2 — s + 1 = 0 of smallest modulus (~ log 2.08102). 
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The iterates of (f ar^ 


Po = [x, y, 

Pi = lzAi,x^y,xyz] 

P2 = [z A2 ,x Al,x^ yAi] 

P:i = [xy ^3, 2;^ Ai x^^y z Aj A2] 

P4 = [x'^y ^1^14,2; A2 A3, X z'^ Al AI A3] 

P 3 = [z A 3 A 2 As, X y"^ A 3 A 4 , y 2 :^Ai Al A 3 A 4 ] 

Pe = Al A3 Ag, xyAiAi A^, x y^ ^A2 A^ A^ Ag] 


The form stabilises into 


Pk — [gi,k Afc _4 Afc _3 Afc, g 2 ^k Afc-s Afc _2 A^_ 3 , g3^k Afc_4 A ^_3 A‘1_2 Afc_i]. (20) 


The coefficients gi^k depend on the initial condition {x,y,z}, and verify gi^k +9 = gi,k- They 
are given by: 


k mod 9 

i = 1 

i = 2 

i = 3 

0 

X 

y 

z 

1 

z 

x^y 

xyz 

2 

z 

X 

9 

X y 

3 

xy 


x^yz 

4 

9 

X y 

z 

9 

xz 

5 

z 

9 

xy 

yz‘^ 

6 

z^ 

xy 

xy‘^z 

7 

y 

z 

9 

xy 

8 

xy"^ 

z 

xyz 


Notice that the presence of the factors gi^k shows that the coordinate planes appear period¬ 
ically in the iterates. 

The recurrence relations between the A’s also have coefficients depending on the initial 
conditions, and on fc in a periodic way (period 9). 


^k ^k-A ^k + Pk ^k-3 ^k-2 ^k-1 + ® ^kPk Ak-A Afc-3 Afc_2 Afc_i Afc — Afc-S A 


k+1 


( 21 ) 


The coefficients are given by the following table 


k mod 9 

A 

p 

0 

X 

z 

1 

1 

xy 

2 

z 

9 

X y 

3 

1 

xz 

4 

xy 

z^ 

5 

1 

yz 

6 


9 

xy 

7 

1 

xy 

8 

y 

z 


^From Ti = BiCi, the blocks Ak are products of two blocks, which we do not write for simplicity 
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The recurrence relation (21) differs from the one given for a = 0 in 30 , since it depends 
on the order k. The explicit calculation of the hrst iterations indicates that it verihes the 
Laurent property for arbitrary x, y, and a. 


3.7 An integrable map in P 3 : N=3 Periodic Volterra 


if := [x, y, z, t] —)■ [x', y\ z\ t'] with 


x' = —x{f‘ + 2et{y — z) — e^x^-\-e^{y + zy‘) (22) 

t' = —t{f‘ + e^{x^ + y‘^ + z^ — 2xy — 2yz — 2xz)) (23) 


y' and z' being obtained from ( p^ by circular permutations of x, y, 

The map comes from a discretisation of a continuous integrable system, and is known to 


have two algebraic invariants 44 


Starting from po = -Bo, Cq, Hq] we get a sequence of points of the form: 


pi = [do di, i?o -Bi, Co Cl, Dq D\\ 

p2 = [do di d 2 ,-Bo-Bi-B 2 , Co Cl C 2 ,-Do ^ 1 -Bi Cl] 

p3 = [do di d 2 da, Do Di D 2 D 3 , Co Cl C 2 C 3 , Do Di D 2 d 2 C 2 ] 

p4 = [do d 2 ds d 4 . Do D 2 D 3 D 4 , Co C 2 C 3 C 4 , Do D 3 da C 3 ] 

p5 = [do da d 4 da, Do D 3 D 4 D 5 , Co C 3 C 4 C 5 , Do Di D 4 d 4 C 4 ] 


The form of the iterates stabilises after three steps, with as slight difference in the structure 
of the last component between odd and even order. 

Denoting the fcth iterate pk as [X^., Y^, Z^, T^.] with 


Xfc — do dA;_2dfc_idfc, Yfc — DoDfc_2Dfc_iDfc, Zk — CoCk-2Ck-iCk 
T2m+1 = DoDid2m,D2mC2m, ^2m+2 = Dod2m+lD2m+lC2m+l, 

the various A, B,C, D^s verify sets of constraints of the form 
Do{Xk + Yk + Zk) — (do + Do + Co) Tfc = 0 

Tk + e(±Xk ±Yk± Zk) - (Do + e(±do ± Dq ± Co)) ak-i f3k 'fk+i = 0 


(24) 


(25) 


where {a,/ 9 , 7 } is some permutation of {d,D,C}. 

These constraints, extending over strings of successive points of length 3, are responsible 
for the factorisation properties: they define various ideals, and the factorisations take place 
in the the algebra generated by the d, D, C,... quotiented by these ideals. Moreover the 
constraints are conserved by the evolution. The various ± signs in (25) depend on fc in a 
periodic way (period 3). 
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A typical example of these relations is (for k = 3): 


Do (A^3 + ^3 + Zo) — {Aq -\- Bq-\- Co) T 3 = 0 
T3 + e(—X3 + 1(3 + Z3) — (Do — eAo + cBq + eCo) A2 B3C1 = 0 
T3 + e(X3 — + Z3) — (Do + eAo — eBo + eCo) Ai B2C3 = 0 

< T3 + e(X3 + Y3 — Z3) — (Do + eAo + eBo — eCo) A3 B1C2 = 0 
T3 + e(—X3 — Y3 + Z3) — (Do — eAo — eBo + eCo) Ai B3C2 = 0 
T3 + e(—X3 + 1^3 — Z3) — (Dq — eAo + eBo — eCo) A3 B2C1 = 0 
T3 + e(X3 — Y3 — Z3) — (Do + eAo — eBo — eCo) A2 B1C3 = 0 


The set of constraints is invariant by circular permutation A ^ B ^ C ^ A. It cannot be 
written solely in term of the components X, X, Z, T. It cannot be solved straightforwardly 
for any set Ak,Bi:,CkS, because it is then over-determined. 

On the other hand, writing that [X^+i, Yfc+i, T^+i] is the image of [X^, X^., T^] 

by ip yields algebraic equations for {A^+i, B^+i, Ck+i}, which one can solve rationally. 

The effect of the set of constraints (25) is that these expressions can be simplihed. The 
factors of {A^+i, Bk+i-, Cfc+i}, reduce to monomials in A^-i, Bk-i, Ck-i,Ak- 2 , Bk- 2 , Ck -2 and 
possibly Di for odd k. 

The simplihed relations dehning the iteration read: 


+ ‘^DTkiYk — Zk) + e^(Yk + Zk)"^ — e^X‘1 + XkA]^j^-iBk-2Ck-2^k-iBk-iCk-i = 0 
TI + 2eTk(Zk — Xfc) -|- e^(Zk + XkY — e^Y^ -|- XkAk-2Q]i+iCk-2^k-iBk-iCk-i = 0 
TI + 2eTk(Xk — Yk) + e^(Xk + X^)^ — e^Z]: -|- AfcAfc_2i?fc-2Ck+iAfc_ii?fc_iC'fc_i = 0 


(26) 


with Afc = Di for odd k and A^ = 1 for even k. 

These equations are linear in Bk+i, Ck+i}- They tell us that the factor 

fk = Afc Ak-2 Bk-2 Ck-2 Afc_i Bk-i Ck-i 

goes away from the homogeneous coordinates when calculating Pk+i as ip(pk)- They extend 
over a string of successive points of length 4. They are the generalisation of the Hirota 
bilinear formalism for the map under consideration, but they are not quadratic anymore. 
They do not have the Laurent property. 


Thanks to the relations (25), their solution in {Ak+i, Bk+i, Ck+i} is polynomial in terms 
of the initial conditions {Ao, Bo, Co, Do}, and it is possible to show that Ak+i, Bk+i,Ck+i 
are the proper transforms of Ak,Bk,Ck- In other words, the A, B, CA do not factorise. 


3.8 A linearisable map 

Linearisable recurrence are known to have special singularity structure. As an example we 


can take the one studied in 45 , where it was shown to be non-conhning, but integrable. 
The recurrence is 


1 

'^n+1 '^n “I" Z ' 

1 Uji Uji—l 
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so that 


if : [x, y, z] — > [2 xz + — xy — yz, x {x — y + z), {x — y + z) z] 

i> :[x,y,z] —y [y {x - y - z) ,xz + xy - 2yz - y"^, z {x - y - z)] 
= z [x — y + z) = z Bl, = z [x — y — z) 

The iterates pk take the form 

Pk = [Ak, Ak-i Bk, z Bi B 2 B^ ... Bk] 

with 


(p*{Bk) = z Bk+i 

There is regularity of the pattern, but the number of factors increases with k. This is 
related to the fact that at each step, one more factor Bi appears, and is at the origin of the 
low (linear) growth of the degrees of the iterates. 

The recurrence relations between blocks read 


z Bi B 2 ... Bk-i (2 Ak Bk — Ak-i Bk — Ak+i) + Ak{Ak — Ak-iBk) — 0 
z Bi B 2 ... Bk-i {Bk — Bk+i) + Ak — Ak-i Bk = 0 


(27) 


The linear growth of the degrees can be read from the previous relations. 


3.9 An unruly model 


We know of maps for which the sequence of degrees does not verify any hnite recurrence 
relation. Although this does not prevent their entropy from being the logarithm of an 
algebraic integer, it will prevent the existence of the Hirota like forms we have seen in the 
previous cases. 

A simple example was found in [^. It is a monomial map in three dimensions: 

p : [x, y, z, t] —^ [ty, tz, x^, tx], t[) ■. [x, y, z, t] —^ [tz, xz, yz, 

= X^t, = Zt^ 


This map has the peculiarity that the entropies of p and of 'ip differ 29 


The structure of the iterates is simple, since they are all written in term of the coordinate 
planes 




for some powers 5*. For ip, the form of the iterates stabilises and sequence of degrees verihes 
a hnite recurrence relation, but this is not the case for p. The peculiarity of the model is 
that the singularity structure is such that the sequences of proper transforms which are at 
the basis of the observation we made for all other examples do not appear here. This model 
is a limiting case to keep in mind for further developments. 
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3.10 A delay-differential / differential difference equation 

Consider the following eqnation: 


a u(t) — b dtu{t) = u{t) {u{t + 1) — u{t — 1)) 


( 28 ) 


where dt means time derivative. 

Equation (28) was obtained in by a non trivial reduction of a semi-discrete equa¬ 
tion This equation is a delay difference equation of which the entropy has been evaluated 
in (48[/and found to be vanishing. One may equivalently consider the differential difference 
equation, or recurrence of order two defined on functional space: 

(X b dfUnit^ Un(t') ('Un+l(^) 1(^)) (29) 


The maps (p and tp associated to these equations are: 

(p-.[x^y^z\ —> [axz — b{x'z — xz')-\-xy^ x'^, xz] (30) 

Ip : [x, y, z] —)■ -ayz + b {y'z - yz') + xy, yz] (31) 


were prime (') means derivative. Here x, y, z should be considered as a container for the 
inhnite sequences [x{t),x'{t),x"{t),...], [y{t),y'{t),y"{t),...], and [z{t), z'{t), z"{t),...]. 

For this map 

= x^, K^{[x,y,z]) =y^. 


and we get the following form for the hrst iterates starting from po 


Po = [^0, Bq, Cq] 

Pi = [ Ai, Hq, AqCq ] 

P2 = [ ^25 Af, AqAiCq ] 

PS = [ ^0^3, A 2 , AqAiA2Cq ] 

Pi = [ AP\Ai, AqA\, AiA2A^Cq ] 

Pk = [ ^fe-3 Ak, Ak-i Ak-2, Ak-2Ak-l Cq ] (32) 


The recurrence for Ak reads 


o Cq Ak Ak+l Ak+2 Ak+3 + Aq Ak+3 A^k+2 + ^ ^‘k^‘k+3 


Cq Ak+1 Ak+2 


Ak A 


k +3 


— ^k+l -^k+4 


(33) 


This relation extends over a string of length 5. Again, although given as a fraction, Ak+i is 
a differential polynomial in the initial conditions. 

The proof of relation (33) is done by recursion. It is verihed for k = d and k = 1. We 


moreover know that the pullback of any Ak by p is of the form AqA^+i, it is easy to show 


the validity of (33) for k + 1. In particular, one hnds that Ak+\ is the proper transform of 
Ak and moreover 


/ Cq Ak+I Afc+2\ 

\ Ak Ak+3 ) 


9a Ak+2 Ak+3 

Ak+1 Ak+i 
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so that the pullback of the derivative term appearing in (33) does not contain factors Aq. 
We also get relations on the various degrees 6{Ak): 

1 + 6{Ak) + (5(y4fc+i) + 6{Akj^2) + = 1 + 2 5(y4fc+2) + <5(y4fc+3) 

= ^{Ak)+ 2 5{Ak+i) + 5{Ak+A)- 

One then easily proves the result on the sequence of degrees of pn given in (48] 


^(Pn) = o 6 + 9 - (-1)’ 


ensuring the vanishing of the entropy. 

One could always question the notion of integrability for differential-difference equations, 
and even more for delay-difference equations, but the vanishing of the algebraic entropy is a 
very strong structural constraint on the equation. 


3.11 An integrable lattice map: Qy 

The model, introduced in |49|[50] is defined on a plane square lattice, by a multilinear relations 
between the values of an unknown function Mn,m, n ^ Z,m ^ Z. It interpolates between the 
various models of the Adler-Bobenko-Suris list 51-53 , and has seven free parameters. The 


integrability of the model was originally based on the evaluation of its algebraic entropy, 
which vanishes. Subsequently this model was shown to have an infinite set of symmetries. 


implemented by two recursion operators related by a elliptic condition 54,55 
The elementary cell of the lattice, written with the usual convention Ur. 

^25 ^n+l,m+l '^12 looks like 


— Un-\-l^Yn — 



and the local relation dehning the model reads: 


CLl U U2UJ2 T ®2 ("W U2 T U U2 U ^2 U Uj V,i 2 T Ux U2 U12) T 03 {u Ui -\- U2 U12') 

-I-CI 4 (it^ U 2 -\- u U 12 ') [u U 2 + 'Uj OiQ {u Uj -\- U 2 -\- U 12 ') + CI 7 = 0 (34) 

The previous relation being multilinear, it is possible to calculate any of the corner variables 
in term of the other three. On each cell set (pn,m ■ u —> Ui 2 and : Ui 2 —> u. 

In order to define an evolution we need to specify initial conditions. We choose to give 
initial conditions on two adjacent diagonals, labelled —1 and 0 , and use the local condition 
to hll the entire lattice. Points on the diagonal k have coordinates {n, m} with n + m = k. 
We may then define a map p from diagonal k to fc -|- 2, and "0 from diagonal A; -|- 2 to fc 
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(straight arrows in Figure 1), for which the values on diagonal fc + 1 enter as parameters. 
Although the space of initial conditions is inhnite, a point {n, m}, ’sees’ only a hnite number 
of initial points on the diagonals —1 and 0. We projectivise the system by turning the space 
of values at each point in to a projective line Pi, with homogeneous coordinates [Xn^m, W,m] 
so that Un,m = Xn,m/yn,m, and Writing only polynomial expressions, keeping in mind that 
any common factor to X and Y ought to be removed. 



Figure 1: Initial conditions and (north-east) evolution 


We know from 50 


that the drop of the degrees of the successive iterates is intimately 

This biquadratic is nothing but the multiplier 


related to one of the biquadratics given in 51 
Kn,m — Calculated ou oue cell: to any pair of adjacent points {n, m} 

on a diagonal we associate the polynomial 


f^n,m {QjiQj'i Clg) Xn^m+lYn,m+lXfi-\-l,mYn-\-l,m 

-|-(aia6 + 0204 ~ ~ 02 ^ 5 ) (-^n,m+l^n+l,m-^n+l,m + W,m+l) 

^ T(c12®7 T U 4 CI 6 ®5®6) (hj^ T (35) 

T(®2®6 ®3®5) {^n,m+lYn+l,m T 

+ (0104 - al ) + (04O7 - al ) 

The key fact is that these polynomials split into two factors as soon as n + m > 2. One 
of these two factors is common to Xn+i,m+i and W+i,m+i the second one is common to 
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Xn+ 2 ,m +2 and Fri+ 2 ,m+ 2 - The evolution equations may be rewritten as the recurrence 


Xn4-l.m4-l = ^2Xn^rn+lXYi^rnXri+l^m 


^A^n^rn-\-l^n-\-l^rri^n^rn ^^^n^rn+l^n^rri^n-\-l^rn ^7^n^rn+^y^n^rn^n+l^rn 

Uq ^XYi^rn+\Yn^rny^n+l^rri T XYi^rny^n^rn-\-\^n-\-l^rn T Xn^-i^^YYi,'m-\-lYn,rn^ 

Yn4-l,m+l ^n+l,m+l ^l^n,m-\-l^n,ni^n-\-l,rn T ^S^n^m+lY^n^rn^n+l^rn 

T ^4^n,rn^n^rn+l^n-\-l,rn T ^5^n+l^rn^n,rn+l^n^rn T ^G^Ujm-l-l^njrn^n-l-l^rn 

-\- CI2 -\- Xfi_^YnXfiJ^i_^rfiYfi^fyf^^i -\- Xfi_^^j^iXfiJ^i_^YYiYfi_^^^ , 

with the added initial condition that i}k,i = 1 for k + I = 1. The various points entering the 
defining relation of the f2’s are pictured in Figure 1 with the curved arrows. The f2’s and 


( 36 ) 


the X, F’s given by (36) are polynomials in the initial conditions, and there is no additional 


common factor to X^+i^^+i o,nd Yn+i^m+i for generic values of the parameters a. Relations 
(36) yield the sequence of degrees of the iterates found in 49,50 , quadratic growth and 


vanishing entropy. More details will be given elsewhere. 


4 Conclusion and perspectives 


We have shown that, for systems undergoing a rational discrete evolution, a self-organisation 
takes place after a hnite number of steps: the structure of the iterates stabilises - this is to 
be compared with the results of - and the r functions pop out spontaneously as pieces 
of the components of the iterates, providing a new set of variables to describe the evolution. 

This change of description of the models is not to be confused with a usual (birational) 
change of coordinates: we barter the original coordinates for pieces of the components of 
strings of successive iterates and transforms of the factors of the multipliers k^. 

The recurrence relations obeyed by the new variables provide us with an exact evaluation 
of the algebraic entropy and of its avatars obtained by reductions to integers and hnite helds, 
without restriction to integrability. They support, by the form they take, the fundamental 
conjecture presented in (T^ that the entropy is always the logarithm of an algebraic integer. 


They also invite us to make contact with the results of 57,58 on orthogonal polynomials. 


In the integrable cases our approach may hnally yield, in addition to the known ap¬ 
plications (Lax pairs, special solutions, soliton solutions, grassmanian description) a new 
classihcation tool. 

All of this is matter for future work. 
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